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Electron transport through a quantum dot coupled to superconducting leads shows a sharp con¬ 
ductance onset when a quantum dot orbital level crosses the superconducting coherence peak of one 
lead. We study superconducting single electron transistors in the weak coupling limit by connect¬ 
ing individual gold nanoparticles with aluminum junctions formed by electromigration. We show 
that the transport features close to the conductance onset threshold can be accurately described 
by the quantum dot levels’ hybridization with the leads, which is strongly enhanced by the diver¬ 
gent density of states at the superconducting gap edge. This highlights the importance of electron 
cotunneling effects in spectroscopies with superconducting probes. 


The transport properties of a tunnel junction to a sin¬ 
gle quantum dot (QD) are a sensitive probe of the its 
spectral properties. As the bias voltage across the junc¬ 
tion is ramped up, new quantized transport channels - 
associated to the QD energy levels - become accessible, 
leading each to a stepwise increase in the current [1-4]. 
The electronic tunneling process associated to each chan¬ 
nel can be elastic - reflecting thereby the QD energy spec¬ 
trum - or inelastic, in which case additional bosonic exci¬ 
tations are involved [5-7] . When normal conductive leads 
are used, thermal broadening limits spectroscopic resolu¬ 
tion to about 3.5in weakly coupled QD junctions. 
The use of superconducting (S) leads strongly modifies 
the picture. As soon as the thermal energy VbT is a 
few times smaller than the superconducting gap A, ther¬ 
mal quasiparticles in the leads have a vanishingly small 
density. Thermal smearing is thereby avoided and the 
spectroscopic resolution can be much improved [7-13]. 
Ultimately, another limitation to spectroscopic resolution 
arises from the tunnel coupling induced hybridization of 
the isolated quantum level at study. With normal con¬ 
tacts, this manifests as a lorentzian broadening of the 
spectral features, set by the tunnel coupling energy scale 
h'y. Nevertheless, this description does not hold in the 
case of superconducting leads [14, 15]. 

In this Letter, we report transport measurements in 
S-QD-S junctions in the weak coupling limit, in a regime 
where a unique separation of all relevant energy scales 
allows for quantitative comparison to theory. We exper¬ 
imentally show that the large DOS at the superconduct¬ 
ing gap edge of the leads strongly enhances hybridiza¬ 
tion of the single quantum level it is tunnel coupled to. 
Hybridization of the QD spectral function is found to 
be the dominating limitation to spectroscopic resolution, 
despite the weak coupling. 

The device fabrication process relies on controlled elec¬ 
tromigration of an on-chip all-metallic device presenting 
a constriction [16]. This technique was successfully ap¬ 
plied for connecting single molecules in numerous exper¬ 
iments [5, 17-21]. Using aluminum constrictions, elec¬ 
tromigration has been used for forming superconducting 


single molecule transistors [22]. In that experiment, the 
use of a narrow section of proximity superconducting gold 
to contact the molecular QD yielded rather high tunnel 
couplings, at the cost however of losing a hard supercon¬ 
ducting gap in the leads. Alternatively, pure aluminum 
leads provide more weakly coupled devices, with BCS 
type superconducting contacts. We pattern aluminum 
constrictions ~ 100 nm wide and 15 nm thick on top of 
a locally defined backgate. Optimal gate coupling along 
with galvanic isolation up to gate voltages ^ 8 V is pro¬ 
vided by a 10 nm thick atomic layer deposition of alumina 
on top the backgate electrode. Gold nanoparticles of 5 
nm diameter [23] are dry-casted onto the sample from a 
toluene suspension shortly before introducing the sample 
into the cryostat. We perform electromigration at 4.2 K 
in the cryogenic vacuum of a home-made inverted dilu¬ 
tion refrigerator. A bias voltage is slowly ramped up 
30 mV/s) while monitoring the constriction conductance. 
A real-time controller sets the bias to zero within mi¬ 
croseconds as soon as the constriction conductance falls 
below 100 fi S. In spite of the high access series resistance 
of 170 U per wire, which is due long lossy coaxial lines at 
base temperature on each lead for filtering, electromigra¬ 
tion typically sets on at about 1 V. Electromigration gaps 
are difficult to see in the Scanning Electron Micrsocope, 
especially in the presence of the spin-coated nanoparti¬ 
cle layer. Figure la shows the image of a large gap with 
virtually infinite tunnel resistance for the purpose of il¬ 
lustration. 

After electromigration, most junctions display tunnel 
behavior with resistances in the upper MU range. Sam¬ 
ples showing stable gate-dependent conductance features 
are further investigated in situ at dilution temperatures. 
All measurements reported here were performed at an 
electronic base temperature T « 100 mK. The current- 
voltage I(V) characteristics are measured in a heavily 
filtered DC transport environment; differential conduc¬ 
tance dl/dV maps at different gate voltages Vq are ob¬ 
tained by numerical differentiation. 

The two samples that we report on here are charac¬ 
terized in Table 1. Both samples have charging energies 
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FIG. 1. (a) Scanning electron micrograph of bare A1 con¬ 

striction after electromigration. The scale bar is 250 nm (no 
nanoparticles present here), (b) Schematic of the S-QD-S de¬ 
vice, introducing the different capacitances and tunnel cou¬ 
plings. (c) dl/dV differential conductance map of the asym¬ 
metric device A, displaying its sole experimentally accessi¬ 
ble degeneracy point. The conductance ranges from about 0 
(black) to 1.6 fiS (bright) (T — 100 mK, zero magnetic field). 
A gap at the degeneracy point is visible, (d) Zoom on the 
degeneracy point at device A at a field of 600 mT. The gap 
has disappeared and the conductance onset is Zeeman split 
at one Coulomb diamond edge. 


U > 100 meV and display the well known transport prop¬ 
erties of weakly coupled S-QD-S transistors [3]. Following 
the orthodox single electron transistor picture [1, 2, 24] 
the device is characterized by two tunnel couplings 7 \ to 
the source and drain leads and three capacitances (C*) 
to the gate (G), source (S) and drain (D) electrodes re¬ 
spectively, see Fig. 1. The dimensionless gate coupling 
parameter a = Cq/Cy, (with Cy = Cq + Cs + Cd ) indi¬ 
cates the efficiency of the conversion of the gate voltage 
applied into variations of the QD chemical potential; val¬ 
ues of a in the 10 -2 range are usual in electromigration 
junctions. Because of the extremely high charging energy 
however, only one charge degeneracy point is accessible 
in each device, at Vq. We name N and N + 1 the num¬ 
ber of electrons on the QD at Vq smaller/larger than Vq 
respectively. All devices we have measured showed an ex¬ 
tra gate independent contribution to the current. We at¬ 
tribute this to a direct superconductor-to-superconductor 
(SIS) transport channel, shunting the QD due to the nar¬ 
rowness of the gap produced by electromigration. This 



Rs 

(Gfi) 

a 

Cd/Cs 

1+ /_ %' y s h'/D 
(pA) (pA) (fieV) (^eV) 

A 

0.5 

0.05 

0.70 

166 

92 5.2 0.4 

S 

1 

0.09 

1.07 

290 

250 2.1 1.4 


TABLE I. Device parameters of devices A and S. R s is the 
shunt resistance, a the gate coupling and Cd/Cs the capac¬ 
itive asymmetry. For the definition of the on-state currents 
(/+ and /_) and the tunnel couplings 7 s,d, see text. 


small contribution corresponding to a shunt resistance 
R s in the Gfl range can be well fitted as a standard SIS 
tunnel junction and subtracted from the data. At the de¬ 
generacy point, a spectroscopic gap of total width 4A « 1 
meV persists in the I(V) traces (Fig. lc and 3a). This 
is because elastic single quasi-particle transport at low 
temperatures through a S-QD-S device requires a min¬ 
ima the filled states of one lead to be above the empty 
states of the opposite lead, that is \V\ > 2 A/e. When a 
magnetic field is applied, the spectroscopic gap gradually 
decreases and disappears above 500 mT (Fig. Id). 

Discrete and sharp peaks in the dl/dV maps demon¬ 
strate that individual single orbital levels are addressed 
when sweeping the bias (Fig. lc). The QD level spac¬ 
ing SE exceeds 1 meV, in agreement with earlier works 
on similar devices with normal leads [25]. At sufficiently 
high magnetic fields, the Zeeman splitting of the QD or¬ 
bital level involved, proportional to the applied magnetic 
field, is well resolved in the conductance data (Fig. Id). 
Because of the strong asymmetry of the tunnel rates to 
both leads of device A, shown in Fig. 1 , this splitting is 
only seen along a single edge of the N occupied Coulomb 
diamond. Since the first excited state of the N N + 1 
transition does show Zeeman splitting, we can conclude 
that N + 1 is odd in this device [3] . The associated elec¬ 
tron Lande factor is 2.4T0.3. This slightly larger value 
than previously reported results for gold nanoparticle 
junctions [25] may be due to the fact that the Zeeman 
energy is determined here at rather small fields. 

In the normal state of the leads, the on-state currents 
of the QD device I± (where db reflects the sign of the ap¬ 
plied bias voltage) are directly related to the tunnel cou¬ 
plings 7 s,d [24]. The same is true with superconducting 
leads at bias voltages well above 2A/e. Notably, both de¬ 
vices presented here have comparable on-state currents, 
although quite different coupling asymmetry (Table 1): 
sample A has quite asymmetric tunnel couplings, as op¬ 
posed to the more symmetric sample S. Further, the tun¬ 
nel couplings, in the fieV range, are slightly smaller than 
ksT ~ 10 fieV. They are also significantly below all other 
energy scales, that is, A « 260 /ieV, which is itself much 
less than the level spacing SE (several meV) and the 
charging energy U > 50 meV. This very well separated 
hierarchy of energy scales is a rather unique situation in 
that it allows for a truly quantitative analysis of charge 
transport through a single quantum level. 
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When coupled to a shapeless continuum, a discrete 
quantum level decomposes into a lorentzian spectral 
function. In the presence of a non-trivial (normalized) 
superconducting density of states p(E) in the leads, given 
by the real part of \E\/y/E 2 — A 2 , hybridization can 
be described via an effective tunnel coupling parame¬ 
ter. Consequently, the tunnel coupling itself depends 
on the relative position of e with respect to the lead’s 
chemical potential, namely 7(e) = fryp(e)/2. A Green’s 
function based calculation yields the hybridized QD spec¬ 
tral function A(E,e) = (1/7 x)Xm [(E — e — E) _1 ], which 
is shown for various bare QD energy levels e in Fig. 
2a. Here Xm represents the imaginary part and E = 
— h^\E\/^/A 2 — E 2 is the self-energy. Marked distortions 
from a lorentzian line shape (of width fry) are manifest as 
e approaches the superconducting gap edge. At E = A, 
the spectral density falls abruptly to zero; bare QD levels 
below the gap edge are virtually unaffected by the tunnel 
coupling. 

This approach allows to describe cotunneling contribu¬ 
tions to charge current which, by truncating higher-order 
Green’s functions and assuming a single spin-degenerate 
level with U +00, is found to be [15] 

4e /* +0 ° 

I{V) = — (2 - n QD ) j [ns(E) - n D {E + eV)] 

(1) 

x _ is(E)jp(E + eV) _ 

(E-e)2 + fi s (E)+j D (E + eV)} 2 

Here the rii are the respective filling factors of the states 
in the source, drain and quantum dot. In the sequential 
tunneling limit (7 \ > 0), the lorentzian functional inside 

the integrand reduces to a Dirac peak and \I(V)\ repli¬ 
cates the leads’ superconducting DOS [3]. Note that we 
have not chosen the usual symmetric biasing convention 
(Vs = —Vd = V/2) but we write Vs = V and Vd = 0. 
This convention has the advantage of being closer to the 
experimental reality. We take into account the explicit 
capacitive dependence of de QD energy eon F, that is 
e(V,V G ) = a(V G -Vg) + VC S /C E . 

The prefactor (2 — uqd) plays an essential role in that 
it accounts for Coulomb interactions on the quantum dot 
[15]. A similar expression to Eq. (1) without this prefac¬ 
tor is found when on-dot interactions are neglected [14]. 
A general expression of the QD filling factor tiqd, ac¬ 
counting for cotunneling and valid for arbitrary bias volt¬ 
age and leads DOS, is given by [26] 


n QD 


2 




7 s ns ~ 7 d n D 
(E - e) 2 + (7s +7£>) 2 



-1 


( 2 ) 

where the functions 7 i and rii in the integrand are eval¬ 
uated at E and E + eV for i = S,D respectively. The 
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FIG. 2. (a) Calculation of the normalized spectral function of 

a single quantum level tunnel coupled (with fry = 8TO -3 A) to 
a superconducting reservoir for three different bare level en¬ 
ergies e (indicated by the position of the circles on the top 
part). The grey line shows p(E) as a reference, (b) Oc¬ 
cupation tiqd of the S-QD-S device as a function of V for 
a(Vb — Vq) — — 65/zeV, as calculated from Eq. (2). We 
assume A = 260 peV and tunnel couplings as for device S 
(Table I). 


spin-degeneracy of the QD level is summed out in the 
above expression. The strong energy dependence of tiqd 
(Fig. 2b) is associated to the leads’ DOS. In the on-state 
(at high bias) the expression of uqd reduces to 


on = o ipnpE-js n s 
QD (1 + n D ) 7£> + (1 + ns)^s ’ 


( 3 ) 


as given by rate equations [24]. Note that in symmetric 
tunnel coupling conditions (7 s = 7 d =7) this expression 
yields n°Q D = 2/3 and, using Eq. (1), I± = ±(2/3)e7. 

The current-voltage characteristics of device S near the 
(TV, TV + 1) degeneracy point are shown in Fig. 3a. As 
expected, \I(V)\ is a smeared replica of the leads super¬ 
conducting density of states. The small asymmetry in 
the capacitive coupling to the leads (see inset Fig. 3a) 
is strongly reflected in amplitude and position of the co¬ 
herence peaks’ replicas. The cotunneling model, leading 
to Eqs. (1) and (2), provides a very good quantitative 
description of the data. In particular, the model cap¬ 
tures well the drain-source asymmetries in terms of I± 
and broadening at the current onset. 

Clear evidence of spectral hybridization is found in the 
behavior of the sub-threshold current measured off charge 
degeneracy. At finite gate detuning, there is a non- 
negligible contribution to current below the threshold 
bias voltage, which is highlighted on a logarithmic scale 
in Fig. 3b. The cotunneling model describes well the 
long tail of the subthreshold current, as well as its sudden 
rapid drop for |W| < 2A/e. The slight theoretical under¬ 
estimation of the subthreshold current can most probably 
be assigned to non-equilibrium effects, both in the leads’ 
and in the QD filling factors [27], that are neglected here. 
The overall shape of the subthreshold cotunneling con¬ 
tribution combined with the current far above threshold 
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FIG. 3. (a) \I(V)\ trace for Vg ~ Vq in the more symmetric 
device S. The line is a fit using the model developed in the 
text. Note the asymmetry of both the on-state currents and 
the shape of the coherence peak replicas at opposite voltages, 
which are well captured by the model. Inset: dl/dV map of 
device S near Vq. Differential conductances are shown red 
when large and positive, white for close to zero and black for 
negative. The dashed lines indicate the gate voltages at which 
the I{V) traces displayed in this figure (a& b) are taken, (b) 
Logarithmic scale representation of |/(V)| trace at Vg — Vq — 
—2.17 mV at the conductance threshold (symbols). Best fits 
assuming cotunneling model (red line), sequential tunneling 
with a 10 fiV gaussian broadening accounting for noise (cyan 
dashed dotted line) and inclusion of a Dynes parameter (77 = 
10 _2 A, dark grey dashed line; 77 = 10 _3 A, light grey dashed 
line) are also shown. 


provide a self-consistent estimate of the tunnel couplings 
to both leads (Table 1), which are the only free parame¬ 
ters. The steep slope of the \I(V) \ curves at V = ±2 A/e 
is independent on the dot level and tunnel coupling and 
reflects the experimental bias noise of ^ 10 /iV. Notably, 
the effect of temperature on the above description only 
appears with weighting ^ exp(— min(A, 5E)/kBT) and 
thermal smearing is thus completely negligible. 

Alternative explanations of the subthreshold current 
behavior could be thought of in terms of noise or an ef¬ 
fective broadening of the leads’ superconducting DOS. A 
sequential-tunneling calculation neglecting hybridization 
and assuming 10 /iV gaussian noise is missing the slow 


lorentzian decay of the current (Fig. 3b). While the in¬ 
trinsic broadening of the aluminum DOS is known to be 
extremely small (< 10 -6 A) [28], we have alternatively 
tried fitting the experimental data assuming a Dynes- 
type correction to the superconducting aluminum den¬ 
sity of states [29] in which an imaginary term ir] is added 
to E in the definition of p(E). The assumption of an un¬ 
realistically large 77 ^ 7 can to some extent account for 
the subthreshold decay of the current (Dynes corrections 
also appear as a lorentzian contribution in this regime). 
A Dynes-type superconducting DOS is however unable 
to account for neither the shape of the coherence peaks 
nor the rapid fall-off at V = ±2A/e (Fig. 3b). This 
confirms that spectral hybridization associated to cotun¬ 
neling is the dominant mechanism at work, accounting 
for all aspects of spectroscopic broadening in our data. 

In conclusion, we have shown that while spectroscopy 
with superconducting probes is known to help avoiding 
thermally induced smearing, the main source of spectro¬ 
scopic broadening is then associated to an enhanced in¬ 
vasiveness of the tunnel coupling, producing hybridiza¬ 
tion. The tunnel coupling to the leads has thus to be 
kept very small for well-resolved spectroscopies, imply¬ 
ing small tunnel currents. This is of particular relevance 
to a panel of cryogenic STM experiments on low dimen¬ 
sional nano structures, in which superconducting tips can 
provide subspectroscopic resolution [7, 12, 13]. 
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